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ABSTRACT: The reflection of a normally incident wideband pulse by a half- 
space whose permittivity and permeability obey the one-resonance Lorentz model 
is calculated. The results are compared with those from frequency-domain reflection 
analysis. In the spectral regime wherein the real parts of the permittivity and the 
permeability are negative, as well as in small adjacent neighborhoods of that regime, 
the time-domain analysis validates the conclusion that it is sufficient to ensure that 
the imaginary part of the refractive index is positive — regardless of the sign of the 
real part of the refractive index. 
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1 Introduction 



Take a half-space occupied by an isotropic, homogeneous, dielectric medium of 
relative permittivity e r (A ) = e^,(A ) +ze"(A ), where A is the free-space wavelength. 
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When a plane wave with electric field phasor of unit amplitude is incident normally 
on this half-space, the amplitude r(A ) of the electric field phasor of the reflected 
plane wave is given by 



where n(Ao) is the refractive index. The foregoing standard result is found in count- 
less many textbooks; see, e.g., Born & Wolf |J and Jackson 0. It can be applied 
to insulators, poor conductors as well as good conductors. 

Let the medium occupying the half-space also have magnetic properties repre- 
sented by the relative permeability /x r (Ao) = /4(Ao) + (Ao)- Then, the previous 
result has to be modified as follows: 



in the context of composite materials that supposedly exhibit {e' r < 0, y! r < 0} and 
{e" > 0, \i!' T > 0} in some spectral regime [[J - the time-dependence exp(— iujt) 
being implicit for all electromagnetic fields, with u as the angular frequency. Shelby 
et al. argued that the real part of n(A ) = n'(A ) + m"(A ) must then be chosen 
to be negative. Consequently, it was concluded |5j that "[pjlane waves appear to 

propagate from plus and minus infinity towards the source [but] clearly 

energy propagates outwards from the source." The conclusion has been employed 
to deduce the amplification of evanescent plane waves, leading to a proposal for 
aberration-free lenses ||. 

Certainly, energy does not flow in the backward direction. Rather, the forward 
direction is defined as the direction of the energy transport, while the phase velocity 




(1) 





This result (|2]) for the reflection coefficient r(Ao) gains particular prominence 
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may then be pointed in the backward direction — a known phenomenon [0. Assum- 
ing the one-resonance Drude model for both e r (A ) and /i r (A ), confining themselves 
to the spectral regime wherein e' r (Xo) < and ^(Xo) < 0, and impedance-matching 
the dielectric-magnetic half-space to free space, Ziolkowski and Heyman |4j] pointed 
out that the reflection coefficient r(Ao) must have finite magnitude. Employing the 
one-resonance Lorentz model for both e r (Ao) and / u r (Ao), McCall et al. |§ recently 
deduced that 

(i) both e^(A ) and yuJ,(A ) need not be negative for 7i'(Aq) to be negative and 
n"(Ao) to be positive, and 

(ii) the magnitude of the reflection coefficient r(Ao) does not then exceed unity. 

Let us emphasize that the spectral regime for n'(Ao) < overspills into the neigh- 
borhoods of the regime wherein both e£.(Ao) and /4(Ao) are negative. 

Whereas frequency-domain analysis for the response of a half-space calls for 
the determination of the signs of n'(Ao) and n"(Ao), the refractive index does not 
enter time-domain analysis. Therefore, time-domain analysis offers an independent 
way of assessing the results of frequency-domain analysis. In this communication, 
we present our conclusions from the response of a dielectric-magnetic half-space 
with negative real refractive index in a certain spectral regime to a normally in- 
cident wideband pulse. The Fourier transform of the reflected pulse confirms the 
frequency-domain calculations of r(Ao) in that spectral regime, and underscores the 
requirement of n"(Ao) > for real materials (which cannot be non-dissipative due 
to causality) to be sufficient. 
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2 Theory in Brief 

Consider the region {z > 0, t > 0} of spatiotemporal space. The half-space z > zg is 
occupied by a homogeneous dielectric-magnetic medium whose constitutive relations 
are given by 



, z > z e , (3) 



D(r,t) = e [E(r,t) + ( Xe *E)(r,t)] 
B(r, t) = fi [H(r, t) + ( X m * H)(r, t)] 

where the asterisk denotes the convolution operation |§ with respect to time, while 
the susceptibility functions 

2ttc . / 2-ncot 

Xe,m(t) = p e ,m *T Bm I — 

x exp f-T7^— ) W(t) (4) 

obey the one-resonance Lorentz model. These susceptibility functions correspond 
to 

= 1 + 1 + (iv /:a e A -V' (5) 



in the frequency domain |TD[. In these expressions, e and /x are the permittivity 



and the permeability of free space; Co = (eoA i o)~ 1//2 is the speed of light in free space; 
U{t) is the unit step function; p e>m denote the so-called oscillator strengths; while 
N e ^ m = 27rM e m and \ e , m determine the resonance wavelengths and the linewidths. 
The region z < Z( is vacuous. 

At time t = 0, an amplitude-modulated wave is supposedly launched normally 
from the plane z = 0; therefore, all fields are independent of x and y. The initial 
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and boundary conditions on E(z, t) and H(z, t) 9X6 9S follows: 
E(^,0) = V^>0 

H(z,0) = Vz>0 

E(0, t) = (/(t) sin y Vt>0 

H(0, t) = -(e //i ) 1/2 sin (go*) x Vt > ^ 

Whereas A car is the carrier wavelength, the function 

c t 



(7) 



g(t) 



2X r 



exp(-c t/A 



(8) 



was chosen to represent the pulse. The cartesian unit vectors are denoted by 
(x, y, z). 

The differential equations to be solved 91T6 9S follows: 
d 



dz 



Ey(z,t) = fJtQ 



—H x {z,t) + —( X m*H x ){z,t) 



—H x (z,t) = e 



d d 
-E y (z,t) + -( Xe *E y )(:.li 



(9) 
(10) 



Their solution was carried out using a finite difference calculus described elsewhere 



in detail |TT| . It suffices to state here that z and t were discretized into segments of 
size Az and At, respectively; derivatives were replaced by central differences, and 
the leapfrog method was used . 
Finally, the Fourier transform 



E y (z, X ) 



E y (z, t) exp(— i— - — t) at 



IF 



't a ^0 

was calculated to determine the spectral contents of the incident and the reflected 
pulses. The parameters z, t a and tb were chosen to capture as much of both pulses as 



possible, following the procedure described elsewhere |TT| . The computed ratio of the 
Fourier transform of the reflected pulse to that of the incident pulse is denoted here 
by ttd(Xo), the subscript TD indicating its emergence from time-domain analysis. 



3 Numerical Results and Discussion 



For the sake of illustration, the following values were selected for the constitutive 
parameters of the dielectric-magnetic half-space: p e = 1, p m = 0.8, A e = 300 nm, 
A m = 320 nm, M e = M m = 100. Thus, eJJ(A ) > and ^(X ) > for all A . 
However, e' r (Ao) is negative for Ao G [212.1, 300] nm, but it is positive for all other 
A ; while yu' r (A ) is negative for A G [238.6, 320] nm, and positive for all other A . 
The definition of the refractive index in (§) suggests two possibilities: Either 

A. n"(X ) is negative for A G [236.1, 316.8] nm and positive elsewhere, consistent 
with the requirement of n'(A ) > VA ; or 

B. n'(Ao) is negative for Ao G [236.1, 316.8] nm and positive elsewhere, consistent 
with the requirement of n"(Ao) > OVAo- 

Thus, our attention had to be focussed on the anomalous spectral regime A G 
[236.1, 316.8] nm. In this regime, the reflection coefficient ?"(Ao) for Possibility A 
is the reciprocal of that for Possibility B. The two possibilities can therefore be 
unambiguously distinguished from one another. 

The carrier wavelength was chosen as X car = 240 nm. The pulse duration is 
3 fs and its 3dB band is Ao G [218, 261] nm. Therefore the anomalous spectral 
regime was substantively covered by our time-domain calculations. The segment 
sizes Az = 5 nm and At = 0.015 fs used by us were adequate for the chosen 
constitutive parameters, but obviously would be totally inadequate in the resonance 
bands of e r (A ) and fi r (X ). 

Possibility A is clearly nonsensical. It implies transport of energy in the half- 
space z > Z£ towards the interface z = Zg. Not surprisingly therefore, @, (||) and 
(|) yielded |r(A )| > 1 for all A G [236.1, 316.8] nm. 



6 



Figure 1 presents the computed values of |r(Ao)| obtained from (0), (Q) and (0) 
for Possibility B (i.e., when n"(\ ) > is guaranteed for all A ). The computed val- 
ues of |?"td(Ao)| are also shown therein. The two sets of magnitudes compare very 
well for Ao < 290 nm. Examination of the refracted pulse also showed that it trans- 
ported energy away from the interface z = Zg, which corroborates the observations 
of Ziolkowski and Heyman [f|. 

Thus, time-domain analysis validates the conclusion that n(Ao) must be selected 
in frequency-domain research in such a way that n"(X ) > — irrespective of the 
sign of n'(Ao). 
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Figure 1: Magnitudes of the reflection coefficient computed using two different methods. 
The frequency-domain equations @, (|5|) and (^|) yield the plot labeled FD, after ensuring 
that n"(X ) is positive for all An- The time-domain analysis, followed by the Fourier 
transformations of the electric field associated with the incident and the reflected pulses, 
as discussed in Section 2, yields the plot labeled TD. The time-domain analysis fails in 
the resonance bands of e r and fi r because the chosen discretization of space and time 
is neccessarily very coarse therein. 
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